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Abstract 

Following recent considerations of a non-zero value for the vacuum energy density 
and the realization that a simple Kantowski-Sachs model might fit the classical tests 
of cosmology, we study the qualitative behavior of three anisotropic and homogeneous 
models: Kantowski-Sachs, Bianchi I and Bianchi III universes, with dust and cosmo- 
logical constant, in order to find out which are physically permitted. 

In fact, these models undergo isotropisation, except for the Kantowski-Sachs model 
(f2fc > 0) with < ^A M and for the Bianchi III (Ofc < 0) with < ^a m , and 
the observations will not be able to distinguish between these models and the standard 
model. 

If we impose that the Universe should be very much isotropic since the last scatter- 
ing epoch [z ~ 1000), meaning that the Universe should have approximately the same 
Hubble parameter in all directions, we are led to a matter density parameter very close 
to the unity at the present time. 
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Lately, the issue of whether or not there is a non-zero value for the vacuum energy 
density, or cosmological constant, has been raised quite often. Even taking the Hubble 
constant to be in the range 60-75 km/s/Mpc it is possible to show U that the standard 
model of flat space with vanishing cosmological constant is ruled out. On the other hand, 
if the classical tests of cosmology are applied to a simple Kantowski-Sachs metric and the 
results compared with those obtained for the standard model, the observations will not be 
able to distinguish between these models if the Hubble parameters along the orthogonal 
directions are assumed to be approximately equal 0. Indeed, as Collins and Hawking || 
point out, the number of cosmological solutions which demonstrate exact isotropy well after 
the big bang origin of the Universe is a small fraction of the set of allowable solutions of 
the cosmological equations. It is therefore prudent to take seriously the possibility that the 
Universe is expanding anisotropically. Note also that in || it is shown that some shear free 
anisotropic models display a FLRW-like behaviour. Taking all this into consideration, we 
discuss the behavior of some homogeneous but anisotropic models with axial symmetry, filled 
with a pressureless perfect fluid (dust) and a non vanishing cosmological constant. For this, 
we will restrict our study to the the metric forms 

fa* = _ c 2 dt 2 + G 2 (t)rfr 2 + fo2(t) ( + v 2 d(j) 2\ (1) 

\1 — kv 2 J 

with the two scale factors a(t) and b(t); k is the curvature index of the 2-dimensional surface 
dv 2 / (1 — kv 2 ) + v 2 dcj) 2 and can take the values +1, 0, —1, giving the following three different 
metrics: Kantowski-Sachs, Bianchi I, and Bianchi III, respectively. 

Einstein equations for the metric ([I]), for which the matter content is in the form of a 
perfect fluid and a cosmological term, A, are then as follows: 

a b b kc 2 n _ . 2 , nS 

2 - I + T2 + -W = 8nG P + Ac ' ( 2 ) 
a o o z o z 

b b 2 kc 2 „» 2 , , 

2 I + 17 + ~OT = -8ttG^ + Ac 2 , 3 

cr cr & 



- + U T + - u T = -8nG^ + Ac 2 , (4) 



a b a b ^p 2 
a o a b c z 
where p is the matter density and p is the (isotropic) pressure of the fluid. Here G is the 
Newton's gravitational constant and c is the speed of light. If we consider a vanishing 
pressure {p = 0), which is compatible with the present conditions for the Universe, the last 
two equations take the form 

n b & kc 2 . 2 

« + M<Lac», (6) 
a b a b 

and Eq.(5) can easily be integrated to give 

6^ _ Mi _ kc^ A 2 

V ~ ~W ~ ~W + 3 ° ' [ ) 

where M\ is a constant of integration. 

The Hubble constants corresponding to the scales a(t) and b(t) are defined by 
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H a = a/a and H b = b/b. 

Using them to introduce the following dimensionless parameters, in analogy wiht which 
it is usually done in the Friedmann-Lemaitre-Robertson- Walker (FLRW) universes, let us 
define 

^#2 = Q M, (8) 



l b 

kc 2 



fife (9) 



and 



Ac 2 

m 2 



= Ov- (10) 



The conservation equation (|7|) can now be rewritten as 

n M -n k + n A = i. (n) 

Now defining the dimensionless variable y = b/bo where bo = b(to), is the angular scale 
factor for the present age of the Universe, and using Eq.flTTD (taken for t = t ), one may 
rewrite Eq.(^) as 

V = ±H bo] jl + fi Mo (i - 1) + n Ao (y 2 - 1), (12) 

where the density parameters fl and H b with the zero subscript, denote as before these 
quantities at the present time to- I n this way, the number of independent parameters have 
been reduced. Substituting Eq.(7) into Eq.(2) gives 



M p - Mi f + §Ac 2 a6 2 
2JM 1 b-kc 2 b 2 + ±cW 



(13) 

s .-=«/.' 

where M p is a constant proportional to the matter in the Universe, 

M p = 8nGpab 2 . (14) 
Using the procedure above, Eq.(13) can be rewritten in the form 

ttp - tt M + 2 n A = 2 ^, (15) 

where 

^From Eq.(2) one may define a matter density parameter 

n = 8nGp 2> (17) 

2H a H b + if b 

just like in FLRW models, and such that Q = 1 when /c = and A = 0, and which is related 
to flp by 



Although Qm is not the matter density parameter, it performs the same important role. We 
emphasize the fact that if for one particular time H a = Hb and ^a = 1, then, by Eqs.(ll), 
(15) and (18), 3Q = Q p = Qm = &>k', and if < ^a <C 1 and Qm = 1, then, Qk = and 

n ~ i. 

Introducing another dimensionless variable x = a/a , Eq.(13) takes the form 



^(l-f)+^A (-l+Xy 2 ) + ^ 



x = H bo 



(19) 



and its number of independent parameters was also reduced, now at the expense of Eq.(15) 
taken for the present time t = to- 

The behavior of y(t) may be carried out looking for the y values where y = 0. This 
analysis was made by Mariano Moles f| for FLRW models, in great detail. 

There are two ^a values which characterize two zones of distinct behavior of scale factor 
b. Starting with condition y = one may obtain 



n 



A 



y 



y 6 -y 



(20) 



If we consider Q^ = fl^ (y), as a function of y, then this function presents a relative 
maximum and a minimum, that we will denote by Q\ c and Q\ M , respectively. The relative 
maximum depends on Qm in the following way: For Qm < 1/2 we have 



n 



3Q 



M 



\ 



Mo 



M 



o 2 



n 



M 



1/3 



1 



A c 



<(n Mo - i)yn 2 Mo - 1 + (1 - n Mo )/n Mo 

for Q Mo > 1/2 the expression is 

'9 + 2n 



+ 



1/3 



3Q Mo cos 



(n 



Mo 



(21) 



(22) 



The relative minimum is done by 



^A A 



-3fi M „ cos 



6 + Air 



i^Mo ~ 1), 



(23) 



where 9 = cos 



-i 



These expressions are limiting zones of the (Q Ao ,Qm ) plane, 

where y = has three or one solutions (for details see || Moles]). The £Ia m expression is 
also defined for O^ > 1/2, but it has the meaning of a maximum only for Qm > 1- The Qa 
less or equal to &a m imposes the recollapse of scale factor b, while greater values produces 
inflexional behaviors for b. The Q\ values greater or equal to Qa c are physicaly "forbidden" 
because they don't reproduce the present Universe (see ||). Obviously, Oa m < ^a c always. 
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Although we are considering anisotropic models, the Eq.(12) is exactly the same as the 
one obtained by || for the homogeneous and isotropic FLRW models. From Eq.(12) and 
Eq.(19) one obtains the differential equation 

dx = ^(l-f)+^Ao(-l + ^ 2 ) + gg 

dy n Mo (l-y) + n Ao (y3-y)+y 

This equation has to comply with the conditions imposed by Eq.(ll) and Eq.(15) eval- 
uated at to- There are some particular values of the parameter for which this equation has 
exact solutions. However, for the majority of the values of the parameter, the solution has 
only been obtained by numerical integration. 

Although we are dealing with anisotropic models, we may admit that at a certain moment 
of time, which we can take as the present time to, the Hubble parameters along the orthogonal 
directions may be assumed to be approximately equal, H a ~ H b . This hypothesis has been 
considered in 0, Henriques] for the case of a Kantowski-Sachs (KS) model. From this study 
was derived the conclusion that the classical tests of cosmology are not at present sufficiently 
accurate to distinguish between a FLRW model and the KS defined in that paper, with 
(H ao ~ Hb ), except for small values of bo. 

The Eq.(24), with H ao = Hf, , has three distinct possible integrations, one for each k 
value, (see Figure 1). There are several density values for which the corresponding curves 




Figure 1: Scale factors relation, that is, the x y dependence for the three models. 

to k — ±1 models have vertical asymptotes on the right of vertical axis, limiting y values 
below. So, we see that for the KS model, the scale factor a(t) starts from infinity if b(t) 
starts from zero. For the Bianchi I model, the scale factors are always proportional or even 
equal. In this situation we don't have an anisotropic model; in fact, we can easily prove that 
this model is isotropic by a properly reparametrization of the coordinates. For the Bianchi 
III model, the scale factor b(t) never starts from zero, but has an initial value different from 
zero when a is null. The following plot shows the zones where each model acts (Figure 2). 

Taking into account the analysis given in 0, we may get the behavior of y(t), since this 
dimensionless parameter obeys the same differential equation for these models and for the 
standard model. Now, going back to Figure 1, one can then determine the x(t) behavior. The 
plotting below summarizes the several possibilities for the three models: Kantowski-Sachs, 
Bianchi I and Bianchi III models, respectively. 
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Figure 2: The Kantowski-Sachs model corresponds to the region above the straight line; the 
Bianchi III model corresponds the region below the straight line; the straight line represents 
the set of region for the Bianchi I model. 




Figure 3: On the left the scale factor for the Kantowski-Sachs model (Qk > 0) when 
f^Ao < ^a m - On the right the scale factor for the Kantowski-Sachs model (fl^ > 0) when 

&A M < ^A < ^A c - 



The present technology allows us to "see" the epoch of last scattering of radiation at a 
redshift of about 1000, i.e., we can actually observe the most distant information that the Uni- 
verse provides. Thus, we observe a great isotropy from the Cosmic Microwave Background 
Radiation (CMBR) because this radiation possesses a near-perfect black body spectrum |J. 
The high level of isotropy from this epoch to our days imposes that the two Hubble factors 
H a and H^ must remain approximately equals from this epoch to the present. In other words, 
we must impose a high isotropy level from the last scattering onwards, in our expressions, 
i.e., 

AH H n - H h 



H a H a 



such that <C 1. We performed several numerical simulations and concluded that the 
sum Qm + ^A must be close to the unity from above for Kantowski-Sachs and from below 
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Figure 5: On the left the scale factor for the Bianchi III model (Qk < 0) when Q\ Q < Q\ A 
On the right the scale factor for the Bianchi III model (Qfc < 0) when Oa.m < ^Ao < ^a c - 



for Bianchi III models^]. We summarize in the table below the AH/H a at z = 1000, for 
Kantowski-Sachs and Bianchi III models, if we supose H ao = Hi 





^M 


^Aq 


tt 




^(%) 


K-S 


1 


7.46 x 10~ 8 


1 - 4.974 x 10~ 8 


+7.46 x lO" 8 


-0.4 x 10" 4 


io- i6 


1 


3.33 x 10- w 


+io- 15 


-1.0 x 10" 4 


B III 


1 - 10~ 8 


9.9 x 10~ 9 


1 - 9.93 x 10~ 9 


-lo- 10 


+9.4 x 10~ 2 


9.99 x 10~ i3 


1 - lO" 1 ^ 


9.997 x 10" 13 


-io- i5 


+ 1.0 x 10" 4 



We excluded the situations in that Q\ is close to unity, because they are not compatible 
with the Q Q values observed today. The two Q values 1 — 6.67 x 10~ 8 and 1 — 9.93 x 10~ 8 for 
Kantowski-Sachs and Bianchi III models, respectively, were chosen so they would be both 
compatible with the estimated value for our Universe matter density and the restriction 
<C 1 on z = 1000. These two values for Q are within the range used in the plots of 



AH 



right hand side figures (3) and (5), respectively. 

■""It is obvious that for Bianchi I model (fW +^A — 1); with our restrictions, we have always AH/H a = 0. 
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Conclusion 



For the Kantowski-Sachs model, we conclude that if the scale factor b(t) starts from 
zero, then the scale factor a(t) will start from infinity and decreases afterwards. When 
Q\ < ^a m , b{t) reaches the maximum value recollapsing after that. So, a(t) will reach a 
relative maximum, when b(t) is maximum, because x(y) has a relative minimum for y < 1 (see 
Figure 1). After that, when b(t) = 0, a(t) goes to infinity again. When J1 Am < fi Ao < ^a c) 
the scale factor b(t) grows indefinitely, giving place to an inflationary scenario. Then, a(t) 
decreases reaching a minimum value, and growing after that indefinitely, and becoming 
proportional to b(t). The initial singularity is of a "cigar" type. 

For the Bianchi I model (£l ko = 0), the scale factors a(t) and b(t) are proportional or 
even equal. Thus, this model turns out to be an isotropic one. However, when f2 Ao < Q\ M , 
a(t) and b(t) reach the maximum and recollapse after that. And when J7 Am < ^a < ^a c > 
a(t) and b(t) grow indefinitely after an inflection. 

For the Bianchi III model (Qk < 0), when Q\ < £l\ M , b(t) starts from an initial non 
vanishing value (b(t = 0) = b > 0), reaching a maximum and recollapsing after that until 
reaches the same value for t — 0. Also, a(t) has a similar behavior, but starts from zero and 
recollapses to zero. When f2 AM < f2 Ao < fl\ c , b(t) starts again from a non vanishing value 
(bo > 0), growing indefinitely with an inflection. In this case, a(t) starts from zero and grows 
indefinitely becoming proportional to b(t). So, the initial singularity is of a "pancake" type. 

In conclusion, these models undergo isotropisation, except for the Kantowski-Sachs model 
(fi fco > 0) with Q Ao < J1 Am and for the Bianchi III (tt ko < 0) with f2 Ao < fi A . If we impose 
that the Universe should be very much isotropic since the last scattering epoch (z ~ 1000), 
meaning that the Universe should have approximately the same Hubble parameter in all 
directions, we are led to a matter density parameter very close to the unity at the present 
time. 
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